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Within this paper, control-relevant vehicle design concepts are examined using a widely 
used 3 DOF (plus flexibility) nonlinear model for the longitudinal dynamics of a generic 
carrot-shaped scramjet powered hypersonic vehicle. Trade studies associated with vehi- 
cle/engine parameters are examined. The impact of parameters on control- relevant static 
properties (e.g. level-flight trimmable region, trim controls, AOA, thrust margin) and dy- 
namic properties (e.g. instability and right half plane zero associated with flight path an- 
gle) are examined. Specific parameters considered include: inlet height, diffuser area ratio, 
lower forebody compression ramp inclination angle, engine location, center of gravity, and 
mass. Vehicle optimizations is also examined. Both static and dynamic considerations are 
addressed. The gap-metric optimized vehicle is obtained to illustrate how this control- 
centric concept can be used to “reduce” scheduling requirements for the final control 
system. A classic inner-outer loop control architecture and methodology is used to shed 
light on how specific vehicle/engine design parameter selections impact control system de- 
sign. In short, the work represents an important first step toward revealing fundamental 
tradeoffs and systematically treating control-relevant vehicle design. 

I. INTRODUCTION AND OVERVIEW 

Motivation. With the historic 2004 scramjet-powered Mach 7 and 10 flights of the X-43A 1-4 , hypersonics 
research has seen a resurgence. This is attributable to the fact that air-breathing hypersonic propulsion is 
viewed as the next critical step toward achieving (1) reliable, affordable, routine access to space, as well as 
(2) global reach vehicles. Both of these objectives have commercial as well as military implications. While 
rocket-based (combined cycle) propulsion systems 5 are needed to reach orbital speeds, they are much more 
expensive to operate because they must carry oxygen. This is particularly costly when traveling at lower 
altitudes through the troposphere (i.e. below 36,152 ft). Current rocket-based systems also do not exhibit 
the desired levels of reliability and flexibility (e.g. airplane like takeoff and landing options). For this reason, 
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much emphasis has been placed on two-stage-to-orbit (TSTO) designs that involve a turbo-ram-scramjet 
combined cycle first stage and a rocket-scramjet second stage. This paper focuses on control challenges 
associated with scramjet-powered hypersonic vehicles. Such vehicles are characterized by significant aero- 
thermo-elastic-propulsion interactions and uncertainty 1-18 . 

Controls- Relevant Hypersonic Vehicle Modeling. The following significant body of work (2005- 
2007) 7-9, 19-28 examines aero-thermo-elastic-propulsion modeling and control issues using a first principles 
nonlinear 3-DOF longitudinal dynamical model which exploits inviscicl compressible oblique shock-expansion 
theory to determine aerodynamic forces and moments, a ID Rayleigh flow scramjet propulsion model with 
a variable geometry inlet, and an Euler-Bernoulli beam based flexible model. The vehicle is 100 ft long with 
weight (density) 6154 lb per foot of depth and has a bending mode at about 21 rad/sec. The controls include: 
elevator, stoichiometrically normalized fuel equivalency ratio (FER), diffuser area ratio (not considered in 
this work), and a canard (not considered in this work). A more complete description of the vehicle model 
can be found in previous works 7,29 . 

More recent modeling efforts have focused on improved propulsion modeling 30,31 that captures precombus- 
tion shocks, dissociation, wall heat transfer, skin friction, fuel-air mixing submodel, and finite-rate chemistry. 
The computational time associated with the enhanced model is significant, thus making it cumbersome for 
control-relevant analysis. The simple ID Rayleigh flow engine model discussed within 7,19,26,29 will be used 
in the current paper. 

Hypersonic Vehicle Control Issues. Within this paper, we exploit the generic carrot-shaped vehicle 
3DOF (plus flexibility) model presented j n 7 - 19 - 26 . 29 . A myriad of issues exist that make control design for 
this hypersonic vehicle a potentially challenging problem: 

• Input/Output Coupling. For this system, velocity control is achieved via the FER input. Flight path 
angle (FPA) control is achieved with the elevator 32 . However, there is significant coupling between 
FER and FPA. 

• Unstable/ Nonminimum Phase. Tail controlled vehicles are characterized by a non-minimum phase 
(right half plane, RHP) zero that is associated with the elevator to FPA map 28 . This RHP zero limits 
the achievable elevator-FPA bandidth. (BW) 33 ~ 35 . In addition, the rearward situated scramjet and eg 
(center of gravity) , implies an inherent pitch-up vehicle instability. This instability requires a minimum 
BW for stabilization 29 . To address these potentially conflicting specifications, one approach has been 
to exploit the addition of a canard 19,32, 36-38 . It is understood, of course, that any canard approach 
would face severe heating, structural, and reliability issues. 

• Varying Dynamic Characteristics. Within 29 , it is shown that the nonlinear model changes significantly 
as a function of the flight condition. Specifically, it is shown that the vehicle pitch-up instability and 
non-minimum phase (NMP) zero vary significantly across the vehicle’s trimmable region. In addition, 
the mass of the vehicle can be varied during a simulation in order to represent fuel consumption. 
Several methods have been presented in the literature to deal with the nonlinear nature of the model. 

Papers addressing modeling issues include: nonlinear modeling of longitudinal dynamics 28 , heating 
effects and flexible dynamics 9,24,39 , FPA dynamics 36 , unsteady and viscous effects 8,20 , and high fidelity 
engine modeling 30,31,40 . 

Papers addressing nonlinear control issues include: control via classic inner-outer loop architecture 41 , 
nonlinear robust /adaptive control 32 , robust linear output feedback 38 , control-oriented modeling 19 , lin- 
ear parameter-varying control of flexible dynamics 42 , saturation prevention 22,43,44 , and thermal choking 
prevention 29, 44 . 

• Uncertain Flexible Modes and Coupling to Propulsion. Flexible dynamics have been captured within 
the model by approximating a free- free Euler-Bernoulli beam using the assumed modes method 24 . 
Three flexible modes are used to approximate the structural dynamics. A damping factor of £ = 0.02 
is assumed. The associated mode frequencies are u>\ = 21.02 rad/sec, u >2 = 50.87 rad/sec, W 3 = 101 
rad/sec. These modes must be adequately addressed within the control system design process. While 
performance can be improved by increasing controller complexity (e.g. higher order notches) 42 , one 
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must be wary of, and careful in dealing with, modal/damping uncertainty issues. This is particularly 
important because structural flexing impacts the bow shock. This, in turn impacts the scramjet’s inlet 
properties, thrust generated, aft body forces, the associated pitching moments, and hence the vehicle’s 
attitude. Given the tight altitude-Mach flight regime - within the air-breathing corridor 5 - that such 
vehicle must operate within, the concern is amplified. In short, one must be careful that the control 
system BW and complexity are properly balanced so that these lightly damped flexible modes are not 
overly excited. 

• Control Saturation Constraints. Control saturation is of particular concern for unstable vehicles such 
as the one under consideration. Two specific saturation nonlinearities are a concern for any control 
system implementation. 

— Maximum Elevator /Canard, Deflection and Instability. FPA is controlled via the elevator/canard 
combination 36 . Because these dynamics are inherently unstable, elevator saturation can result 
in instability 43 . Classical anti- windup methods may be inadequate to address the associated is- 
sues - particularly when the vehicle is open loop unstable. The constraint enforcement method 
within 43,45 and generalized predictive control 46 have been used to address such issues. 

It should be noted that control surface/ actuator rate limits must also be properly addressed by 
the control system in order to avoid instability. 

— Thermal Choking/Unity FER: State Dependent Constraint. As heat is added within the combus- 
tor, the supersonic air flow is slowed. If enough heat is added, the combustor exit Mach number 
will approach unity, and the flow is said to be thermally choked 47 . If additional heat is added, 
the upstream conditions can be altered. This can (in principle) lead to engine unstart 5 - a highly 
undesirable condition. The amount of FER that causes thermal choking at a particular flight 
condition is referred to as the thermal choking FER, or FERtc- In general, FERtc depends 
upon the free-stream Mach, free-stream temperature, pressure, and density (which depend on the 
altitude), and the flow turn angle (vehicle geometry + AO A + elastic deflection) 29,46 . In addition, 
since the model does not capture what happens when FER > l 28 , it is natural to restrict FER 
below unity. Given the above, it follows that the minimum of these two constraints dictates the 
available FER at a given flight condition. The resulting state dependent FER constraint can be 
computed (on-line) based on the flight condition, and must be accounted for by the control law. 
Here, uncertainty is of great concern because of the potential unstart issues - issues not captured 
within the model. Engineers, of course, would try to “build-in protection” so that this is avoided. 
As such, engineers are forced to tradeoff operational envelop for enhance unstart protection. 

Control- Relevant Vehicle Design Issues. Despite the successful integrated approach taken by the X- 
43A team, as well as other prior successful flight control efforts, far too often aerospace vehicle design has not 
significantly involved the discipline of controls until very late in the vehicle design process or even afterwards. 
Research programs over the past two decades have suggested that for the anticipated hypersonic vehicles, the 
traditional “sequential” approach is not likely to work. This is attributable, in part, to complex uncertain 
nonlinear coupled unstable, non-minimum phase, flexible dynamics together with stringent flight corridor 
and variable constraints (e.g. specific impulse, fuel use, maximum dynamic pressure, engine temperatures 
and pressures). For such vehicles, an integrated multidisciplinary “parallel” approach - involving multiple 
disciplines up front - is essential. This is particularly true when tight flight control specifications must be 
satisfied in the presence of significant uncertainty. 

Goals and Contributions of Paper. This paper addresses a myriad of issues that are of concern to both 
vehicle and control system designers. In short, this paper represents a step toward answering the following 
critical control-relevant vehicle design questions: 

1. How do vehicle/engine design properties impact a vehicle’s static and dynamic properties? 

2. How do these impact control system design? 

3. How should a hypersonic vehicle be designed to permit/facilitate the development of an adequately 
robust control system? 

4. What fundamental tradeoffs exist between vehicle design objectives and vehicle control objectives? 
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More specifically, in this paper we consider how the following parameters impact the static and dynamic 
properties of a vehicle: 

• engine inlet height, diffuser area ratio, compression ramp inclination, engine location (distance behind 
vehicle nose), vehicle eg (center-of-gravity) , and vehicle mass 

Vehicle optimization is also considered. It is specifically shown that a gap-optimized vehicle can “reduce” 
control system scheduling requirements. A classic decentralized inner-outer loop control system architecture 
is used to illustrate how vehicle/engine parameter selection The gap metric represents a system-theoretic 
measure that quantifies the “distance” between two dynamical systems and whether or not a common con- 
troller can be deployed for the systems under consideration 48,49 . Within this paper, the gap metric is used to 
obtain a “gap-optimized vehicle” which “reduces” how much the vehicle varies throughout the trimmable re- 
gion is obtained. A nonlinear pull-up maneuver is used to show that a “gap-optimized vehicle” can “reduce” 
control system scheduling requirements. Future work will examine the utility of pursuing gap-optimized 
vehicles or optimizing vehicles subject to gap constraints. 

In short, this paper illustrates fundamental tradeoffs that vehicle and control system designers should jointly 
consider during the early stages of vehicle conceptualization/design. The paper also sheds light on how 
specific vehicle/engine parameter selections impact control system design - thus providing a contribution to 
control-relevant vehicle design. While vehicle designers may want to use a higher fidelity model (e.g. Euler 
based CFD with boundary layer reconstruction or Navier-Stokes based CFD 50 ) to conduct more accurate 
vehicle trade studies, this paper shows that a (first principles) 3DOF nonlinear engineering model - such 
as that used in the paper - may be very useful during the early stages of vehicle conceptualization and design. 

Organization of Paper. The remainder of the paper is organized as follows. 

• Section II provides an overview of the dynamical model to be used in our studies. 

• Section III presents engine parameter trade study results as well as a new set of nominal engine 
parameter values. 

• Section IV presents vehicle parameter trade study results. 

• Section V presents vehicle optimization results. 

• Section VI discusses how control system design is impacted by vehicle/engine design parameter selec- 
tion. 

• Section VII summarizes the paper and presents directions for future research. 

II. DESCRIPTION OF NONLINEAR MODEL 

In this paper, we consider a first principles nonlinear 3-DOF dynamical model for the longitudinal dy- 
namics of a generic scramjet-powered hypersonic vehicle 7 ~ 9, 19-28 . The vehicle is 100 ft long with weight 
(density) 6,154 lb per foot of depth and has a bending mode at about 22 rad/sec. The controls include: 
elevator, stoichiometrically normalized fuel equivalency ratio (FER), diffuser area ratio (not considered in 
our work), and a canard. The vehicle may be visualized as shown in Figure l 8 . 

Modeling Approach. The following summarizes the modeling approach that has been used. 

• Aerodynamics. Pressure distributions are computed using inviscid compressible oblique-shock and 
Prandtl-Meyer expansion theory 10, 16, 28, 47 . Air is assumed to be calorically perfect; i.e. constant specific 
heats and specific heat ratio 7 = f ^ = 1.4 10,47 . A standard atmosphere is used. 

Viscous drag effects (i.e. an analytical skin friction model) are captured using Eckerts temperature 
reference method 8,10 . This relies on using the incompressible turbulent skin friction coefficient formula 
for a flat plate at a reference temperature. Of central importance to this method is the so-called wall 
temperature used. The model assumes a nominal wall temperature of 2500°I? 8 . While our analysis has 
shown that this assumption is reasonable for conducting preliminary trade studies, the wall temperature 
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used should (in general) depend upon the flight condition being examined. As such, modeling heat 
transfer to the vehicle via parabolic heat equation partial differential equations (pdes) as well as 
modeling a suitable thermal protection system is essential for obtaining insight into wall temperature 
selection 9 . This will be addressed more comprehensively in a subsequent publication. 

Unsteady effects (e.g. due to rotation and flexing) are captured using linear piston theory 8,51 . The idea 
here is that flow velocities induce pressures just as the pressure exerted by a piston on a fluid induces 
a velocity. 

• Propulsion. A single (long) forebody compression ramp provides conditions to the rear-shifted scramjet 
inlet. The inlet is a variable geometry inlet (variable geometry is not exploited in our work). 

The model assumes the presence of an (infinitely fast) cowl door which uses AOA to achieve shock- 
on-lip conditions (assuming no forebody flexing). Forebody flexing, however, results in air mass flow 
spillage 28 . At the design cruise condition, the bow shock impinges on the engine inlet (assuming no 
flexing). At speeds below the design-flight condition and/or larger flow turning angles, the cowl moves 
forward to capture the shock. At larger speeds and/or smaller flow turning angles, the bow shock is 
swallowed by the engine. In either case, there is a shock reflected from the cowl or within the inlet 
(i.e. we have a bow shock reflection). This reflected shock further slows down the flow and steers it 
into the engine. It should be noted that shock-shock interactions are not modeled. For example, at 
larger speeds and smaller flow turning angles there is a shock off of the inlet lip. This shock interacts 
with the bow shock. This interaction is not captured in the model. 

The model uses liquid hydrogen (LH2) as the fuel. It is assumed that fuel mass flow is negligible 
compared to the air mass flow. The model also captures linear fuel depletion. Thrust is linearly 
related to FER for all expected FER values. For large FER values, the thrust levels off. In practice, 
when FER > 1, the result is decreased thrust. This phenomena 28 is not captured in the model. As 
such, control designs based on this nonlinear model (or derived linear models) should try to maintain 
FER below unity. 

The model also captures thermal choking. In what follows, we show how to compute the FER required 
to induce thermal choking as well as the so-called thermal choking FER margin. The above will lead to a 
useful FER margin definition - one that is useful for the design of control systems for scramjet-powered 
hypersonic vehicles. 

Finally, it should be noted that the model offers the capability for addressing linear fuel depletion. 
This feature was exploited for the nonlinear simulation presented in this paper. 

• Structural. A single free-free Euler-Bernoulli beam partial differential equation (infinite dimensional 
pde) model is used to capture vehicle elasticity. As such, out-of-plane loading, torsion, and Timoshenko 
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effects are neglected. The assumed modes method (based on a global basis) is used to obtain natural 
frequencies, mode shapes, and finite-dimensional approximants. This results in a model whereby the 
rigid body dynamics influence the flexible dynamics through generalized forces. This is in contrast to 
the model described within [28] which uses fore and aft cantilever beams (clamped at the center of 
gravity) and leads to the rigid body modes being inertially coupled to the flexible modes (i.e. rigid 
body modes directly excite flexible modes). Within the current model, forebody deflections influence 
the rigid body dynamics via the bow shock which influences engine inlet conditions, thrust, lift, drag, 
and moment 24 . Aftbody deflections influence the AOA seen by the elevator. As such, flexible modes 
influence the rigid body dynamics. 

The nominal vehicle is 100 ft long. The associated beam model is assumed to be made of titanium. 
It is 100 ft long, 9.6 inches high, and 1 ft wide (deep). This results in the nominal modal frequencies 
<jj\ = 21.02 rad/sec, oq = 50.87 rad/sec, W3 = 101 rad/sec. When the height is reduced to 6 inches, 
then we obtain the following reduced modal frequencies: uq = 10.38 rad/sec, uq = 25.13 rad/sec, 
o>3 = 49.89 rad/sec. Future work will examine vehicle mass-flexibility-control trade studies 19 . 

• Actuator Dynamics. Simple first order actuator models (contained within the original model) were used 
in each of the control channels: elevator- jpkj,FER- can ard - (Note: canard not used in our 
study). These dynamics did not prove to be critical in our study. An elevator saturation of ±30° was 
used. 22,43 It should be noted, however, that these limits were never reached in our studies 41 . Within 
this paper, we consider a pull up maneuver that does not result in elevator saturation. Future work 
will consider more aggressive pull up maneuvers where elevator position and rate saturation become 
very important given the vehicle’s (open loop) unstable dynamics. A (state dependent) saturation 
level - associated with FER (e.g. thermal choking and unity FER) - was also directly addressed 41 . This 
(velocity bandwidth limiting) nonlinearity is discussed below. 

Generally speaking, the vehicle exhibits unstable non-minimum phase dynamics with nonlinear aero-elastic- 
propulsion coupling and critical (state dependent) FER constraints. The model contains 11 states: 5 rigid 
body states (speed, pitch, pitch rate, AOA, altitude) and 6 flexible states. 

Unmodeled Phenomena/Effects. All models possess fundamental limitations. Realizing model limita- 
tions is crucial in order to avoid model misuse. Given this, we now provide a (somewhat lengthy) list of 
phenomena/effects that are not captured within the above nonlinear model. (For reference purposes, flow 
physics effects and modeling requirements for the X-43A are summarized within [52].) 

• Dynamics. The above model does not capture longitudinal-lateral coupling and dynamics 53 and the 
associated 6DOF effects. 

• Aerodynamics. Aerodynamic phenomena/effects not captured in the model include the following: 
boundary layer growth, displacement thickness, viscous interaction, entropy and vorticity effects, lam- 
inar versus turbulent flow, flow separation, high temperature and real gas effects (e.g. caloric im- 
perfection, electronic excitation, thermal imperfection, chemical reactions such as O2 dissociation) 10 , 
non-standard atmosphere (e.g. troposphere, stratosphere), unsteady atmospheric effects 6 , 3D effects, 
aerodynamic load limits. 

• Propulsion. Propulsion phenomena/effects not captured in the model include the following: cowl 
door dynamics, multiple forebody compression ramps (e.g. three on X-43A 54,55 ), forebody boundary 
layer transition and turbulent flow to inlet 54,55 , diffuser losses, shock interactions, internal shock effects, 
diffuser-combustor interactions, fuel injection and mixing, flame holding, engine ignition via pyrophoric 
silane 3 (requires finite-rate chemistry; cannot be predicted via equilibrium methods 56 , finite-rate chem- 
istry and the associated thrust- AO A-Mach-FER sensitivity effects 31 , internal and external nozzle losses, 
thermal choking induced phenomena (2D and 3D) and unstart, exhaust plume characteristics, cowl 
door dynamics, combined cycle issues 5 . 

Within [31], a higher fidelity propulsion model is presented which addresses internal shock effects, 
diffuser-combustor interaction, finite-rate chemistry and the associated thrust- AO A-Mach-FER sensi- 
tivity effects. While the nominal Rayleigh-based model (considered here) exhibits increasing thrust- 
AOA sensitivity with increasing AOA, the more complex model in 31 exhibits reduced thrust-AOA 
sensitivity with increasing AOA - a behavior attributed to finite-chemistry effects. 
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Future work will examine the impact of internal engine losses, high temperature gas effects, and 
nozzle/plume issues. 

• Structures. Structural phenomena/effects not captured in the model include the following: out of plane 
and torsional effects, internal structural layout, unsteady thermo-elastic heating effects, aerodynamic 
heating due to shock impingement, distinct material properties, 57 and aero-servo-elasticity 58,59 . 

— Heating- Flexibility Issues. Finally, it should be noted that Bolender and Doman have addressed a 
variety of effects in their publications. For example, within [9,24] the authors address the impact 
of heating on (longitudinal) structural mode frequencies and mode shapes. 

Within [9], the authors consider a sustained two hour straight and level cruise at Mach 8, 85 kft. 
It is assumed that no fuel is consumed (to focus on the impact of heat addition). The paper 
assumes the presence of a thermal protection system (TPS) consisting of a PM2000 honeycomb 
outer skin followed by a layer of silicon dioxide (Si02) insulation. The vehicle - modeled by a 
titanium beam - is assumed to be insulated from the cryogenic fuel. The heat rate is computed via 
classic heat transfer equations that depend on speed (Mach), altitude (density), and the thermal 
properties of the TPS materials as well as air - convection and radiation at the air-PM2000 surface, 
conduction within the three TPS materials. The initial temperature of all three TPS materials 
was set to 559. 67°!? = 100°F). The maximum heat rate (achieved at the flight’s inception) 
was approximately 12 (1 foot aft of the nose). By the end of the two hour level flight, 
the average temperature within the titanium increased by 125 °R and it was observed that the 
vehicle’s (longitudinal) structural frequencies did not change appreciably (< 2%) [9, page 18]. 
When one assumes a constant 15jp^ heat rate at the air-PM2000 surface (same initial TPS 
temperature of 559.67°!? = 100°F), then after two hours of level flight the average temperature 
within the titanium increased by 200°!? [9, page 19]. In such a case, it can be shown that the 
vehicle’s (longitudinal) structural frequencies do not change appreciably (< 3%). This high heat 
rate scenario gives one an idea by how much the flexible mode frequencies can change by. Such 
information is critical in order to suitably adapt/schedule the flight control system. 
Comprehensive heating-mass-flexibility-control studies will be examined further in a subsequent 
publication. 

• Actuator Dynamics. Future work will examine the impact of actuators that are rate limited; e.g. ele- 
vator, fuel pump. 

It should be emphasized that the above list is only a partial list. If one needs fidelity at high Mach numbers, 
then many other phenomena become important; e.g. O 2 dissociation 10 . 


Longitudinal Dynamics. The equations of motion for the 3DOF flexible vehicle are given as follows: 
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where L denotes lift, T denotes engine thrust, D denotes drag, M. is the pitching moment, TV, denotes 
generalized forces, £ demotes flexible mode damping factor, 07 denotes flexible mode undamped natural 
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frequencies, m denotes the vehicle’s total mass, I yy is the pitch axis moment of inertia, go is the acceleration 
due to gravity at sea level, and Re is the radius of the Earth. 

• States. Vehicle states include: velocity v, FPA 7 , altitude h, pitch rate g, pitch angle 6 , and the flexible 
body states ^ 71 , 771 , 772 , y/ 2 , 773 , 773 . These eleven (11) states are summarized in Table 1. 


# 

Symbol 

Description 

Units 

1 

V 

speed 

kft/sec 

2 

7 

flight path angle 

deg 

3 

a 

angle-of-attack (AOA) 

deg 

4 

q 

pitch rate 

deg/sec 

5 

h 

altitude 

ft 

6 

m 

1 st flex mode 

- 

7 

m 

1 st flex mode rate 

- 

8 

72 

2 nd flex mode 

- 

9 

m 

2 nd flex mode rate 

- 

10 

73 

3 rd flex mode 

- 

11 

73 

3 rd flex mode rate 

- 


Table 1. States for Hypersonic Vehicle Model 


• Controls. The vehicle has three (3) control inputs: a rearward situated elevator <5 e , a forward situated 
canard S c a , and stoichiometrically normalized fuel equivalence ratio (FER). These control inputs are 
summarized in Table 2. In this paper, we will only consider elevator and FER; i.e. the canard has been 
removed. 


# 

Symbol 

Description 

Units 

1 

FER 

stoichiometrically normalized fuel equivalence ratio 

- 

2 

5e 

elevator deflection 

deg 

3 

5c 

canard deflection 

deg 


Table 2. Controls for Hypersonic Vehicle Model 


In the above model, we note that the rigid body motion impacts the flexible dynamics through the generalized 
forces. As discussed earlier, the flexible dynamics impact the rigid body motion through thrust, lift, drag, 
and moment. Nominal model parameter values for the vehicle under consideration are given in Table 3. 
Additional details about the model may be found within the following references '~ 9, 19-28 . 

Scramjet Model. The scramjet engine model is that used in 28,60 . It consists of an inlet, an isentropic 
diffuser, a ID Rayleigh flow combustor (frictionless duct with heat addition 47 ), and an isentropic internal 
nozzle. A single (long) forebody compression ramp provides conditions to the rear-shifted scramjet inlet. 
Although the model supports a variable geometry inlet, we will not be exploiting variable geometry in this 
paper; i.e. diffuser area ratio Ad = f ^ will be fixed with Ad = 1, see Figure 2). 

Bow Shock Conditions. A bow shock will occur provided that the flow deflection angle 5 S is positive; i.e. 

5 S = f AOA + forebody flexing angle + Tu > 0° (9) 

where Tu = 6.2° is the lower forebody wedge angle (see Figure 1). If S s < 0, a Prandtl-Meyer expansion 
will occur. Given the above, a bow shock occurs when the following flow turning angle (FTA) condition is 
satisfied: 

FTA = AOA + forebody flexing angle > —6.2°. (10) 

a In this paper, we have removed the canard. Future work will examine the potential utility of a canard as well as its viability. 
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Parameter 

Nominal Value 

Parameter 

Nominal Value 

Total Length (L) 

100 ft 

Lower forebody angle (til) 

6.2° 

Forebody Length (Li) 

47 ft 

Tail angle (r 2 ) 

14.342° 

Aftbody Length (L 2 ) 

33 ft 

Mass per unit width 

191.3024 slugs/ft 

Engine Length 

20 ft 

Weight per unit width 

6,154.1 lbs/ft 

Engine inlet height li., 

3.25 ft 

Mean Elasticity Modulus 

8.6482 x 10 7 psi 

Upper forebody angle (th/) 

3° 

Moment of Inertia l yy 

86,723 slugs ft 2 /ft 

Elevator position 

(-85, -3. 5) ft 

Center of gravity 

(-55,0) ft 

Diffuser exit/inlet area ratio 

1 

Elevator Area 

17 ft 2 

Titanium Thickness 

9.6 in 

Nozzle exit/inlet area ratio 

6.35 

First Flex. Mode (w ni ) 

22.2 rad/s 

Second Flex. Mode (w n2 ) 

48.1 rad/s 

Third Flex. Mode (w n3 ) 

94.8 rad/s 

Flex. Mode Damping (£) 

0.02 


Table 3. Vehicle Nominal Parameter Values 



Figure 2. Schematic of Scramjet Engine 


Properties Across Bow Shock. Let ( M <*,, T^, p^) denote the free-stream Mach, temperature, and 
pressure. Let 7 = f ^ = 1.4 denote the specific heat ratio for air - assumed constant in the model; i.e. air is 
calorically perfect. 9 10 The shock wave angle 9 S = 9 s (M aai S s ,j) can be found as the middle root (weak shock 
solution) of the following shock angle polynomial 28,47 : 

sin 6 0 s + bsin 4 9 s + csin 2 9 s + d = 0 (11) 


where 


b 


Ml + 2 

Ml 


jsin 2 6 s 


2 Ml + 1 

Ml 


' (7+ l) 2 
4 


7 - 1 ' 

Ml _ 


sin 2 S s 



(12) 


The above can be addressed by solving the associated cubic in sin 2 9 s . A direct solution is possible if 
Emanuel’s 2001 method is used [47, page 143]. 


After determining the shock wave angle 9 Sl one can determine properties across the bow shock using classic 
relations from compressible flow [47, page 135]; i.e. M s ,T s ,p s - functions of ( , <5 S , 7) : 


Too 

Ps 

Poo 


M 2 sin 2 (9 s — S s ) 


(27 Ml sin 2 d s + l- 7) ((7 - 1 )Ml sin 2 9 S + 2) 
(7 + l) 2 Ml sin 2 9 S 

1 + (Ml sin 2 9 S - 1) 

7+1 

Ml sin 2 fl s ( 7 - 1) + 2 
27 Ml sin 2 9 S — (7 - 1) 


(13) 

(14) 

(15) 
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It should be noted that for large M^, the computed temperature T s across the shock will be larger than it 
should be because our assumption that air is calorically perfect (i.e. constant specific heats) does not capture 
other forms of energy absorption; e.g. electronic excitation and chemical reactions 10 . 

Properties Across Prandtl-Meyer Expansion. An expansion fan occurs when there is a flow over a 
convex corner; i.e. flow turns away from itself. More specifically to the bow, if S s < 0 a Prandtl-Meyer 
expansion will occur. To determine the properties across the expansion, let (Moo, Too, Poo) denote the 
free-stream (supersonic) Mach, temperature, and pressure, respectively. If we let S = — S s > 0 denote the 
expansion ramp angle (in radians), the properties across the expansion fan (M e , T e , p e ) can be calculated 
as follows 28,47 : 


/7 + 1 -1 

v 1 = a tan 

V 7- 1 

^2 = vi + 5 

f(M e ) = ^ 

h +1 tan ” 1 
/ 7 - 1 

Pe_ 

'l + ^i Ml 

Poo 

.l+V^e 2 . 

Te_ 


Too 

.l+V^e 2 . 




— tan 1 


T) 


(\Zt +T ( M e2 ~ 1 ' ) ) tan 1 (y M e - i) -^2=0 

7 

7-1 


(16) 

(17) 

(18) 

(19) 

(20) 


V\ is the angle for which a Mach 1 flow must be expanded to attain the free stream Mach. 

Translating Cowl Door. The model assumes the presence of an (infinitely fast) translating cowl door which 
uses AOA to achieve shock-on-lip conditions (assuming no forebody flexing). Forebody flexing, however, 
results in an oscillatory bow shock and air mass flow spillage 28 . A bow shock reflection (off of the cowl or 
inside the inlet) further slows down the flow and steers it into the engine. Shock-shock interactions are not 
modeled. 

• Impact of Having No Cowl Door. Associated with a translating cowl door are potentially very 
severe heating issues. For our vehicle, the translating cowl door can extend a great deal. For example, 
at Mach 5.5, 70kft, the trim FTA is 1.8° and the cowl door extends 14.1 ft. Of particular concern, due 
to practical cowl door heating/structural issues, is what happens when the cowl door is over extended 
through the bow shock. This occurs, for example, when structural flexing results in a smaller FTA (and 
hence a smaller bow shock angle) than assumed by the rigid-body shock-on-lip cowl door extension 
calculation. This is certainly a major concern. It leads one to ask the question: What happens to 
the vehicle properties if no cowl door is present? When the FTA is large or when the vehicle Mach is 
low, the shock angle increases and more air mass spillage would occur. Our analysis shows that the 
impact of neglecting the cowl door on the vehicle’s static properties is significant while the impact on 
the vehicle’s dynamic properties is negligible. This will receive further examination in a subsequent 
publication. 

Inlet Properties. The bow reflection turns the flow parallel into the scramjet engine 28 . The oblique shock 
relations are implemented again, using M s as the free-stream input, 5\ — t\i as the flow deflection angle 
to obtain the shock angle 8 1 = and the inlet (or diffuser entrance) properties: Mi, Ti, pi - 

functions of (M a , 8\, 7 ). 

Diffuser Exit-Combustor Entrance Properties. The diffuser is assumed to be isentropic. The combus- 
tor entrance properties are therefore found using the formulae in 28 , [47, pp. 103-104] - M2 = M2 (Mi, Ad, 7 ), 
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T-2 = T 2 (Mi,M 2 ,7), p 2 = P2(M 1 ,M 2 ,'y): 

(21) 

( 22 ) 

(23) 

where Ad = f is the diffuser area ratio. Also, one can determine the total temperature T <2 = T t2 (T 2 ,M 2 , 7 ) 
at the combustor entrance can be found using [47, page 80]: 

T t 2 = 1 +^M| T 2 . (24) 

Since Ad = 1 in the model, it follows that M 2 = Mi, T 2 = T] , p 2 = pi, and T t2 = [l + A^M 2 ] T\ = T tl . 

FER. The model uses liquid hydrogen (LH2) as the fuel. If / denotes fuel-to-air ratio and / st denotes 
stoichiometric fuel-to-air ratio, then the stoichiometrically normalized fuel equivalency ratio is given by 
FER = f y“q 5 - 28 FER is the engine control. While FER is primarily associated with the vehicle velocity, its 
impact on FPA is significant (since engine is situated below vehicle eg). This coupling will receive further 
examination in what follows. 

Combustor Exit Properties. In this model, we have a constant area combustor where the combustion 
process is captured via heat addition. To determine the combustor exit properties, one first determines the 
change in total temperature across the combustor 28 : 

AT c = AT c (T t2 ,FER,H f ,r) c , Cp J st )= x - ^ 2 ) (25) 

where Hf = 51, 500 BTU/lbm is the heat of reaction for liquid hydrogen (LH2), rj c = 0.9 is the combustion 
efficiency, c p = 0.24 BTU/lbm°R is the specific heat of air at constant pressure, and f st = 0.0291 is the 
stoichiometric fuel-to-air ratio for LH2 5 . Given the above, the Mach M 3 , temperature T 3 , and pressure p 3 at 
the combustor exit are determined by the following classic ID Rayleigh flow relationships 28 , [47, pp. 103-104]: 

mi [ 1 + j(7 - mi\ 

(7M3 2 + 1)2 

T 3 = 

P3 = 

Given the above, one can then try to solve equation (26) for A /3 = A /3 ^A/ 2 , 7 ^ ■ This will have a 

solution provided that M 2 is not too small, AT C is not too large (i.e. FER is not too large or T 2 is not too 
small. See discussion below. 

Thermal Choking FER (M 3 = 1). Once the change in total temperature AT C = A T c (T t2 , FER, Flj, r] c , c p , f s t) 
across the combustor has been computed, it can be substituted into equation (26) and one can “try” to solve 
for A/ 3 . Since the left hand side of equation (26) lies between 0 (for A /3 = 0) and 0.2083 (for A /3 = 1), it 
follows that if the right hand side of equation (26) is above 0.2083 then no solution for A /3 exists. Since 
the first term on the right hand side of equation (26) also lies between 0 and 0.2083, it follows that this 
occurs when AT C is too large; i.e. too much heat is added into the combustor or too high an FER. In short, 
a solution A /3 will exist provided that FER is not too large, T 2 is not too small (i.e. altitude not too high), 
and the combustor entrance Mach M 2 is not too small (i.e. FTA not too large). When A /3 = 1, a condition 
referred to as thermal choking 5, 47 is said to exist. The FER that produces this we call the thermal choking 
FER - denoted FERtc- In general, FERtc w iU be a function of the following: M 0 0 , T a 0 , and FTA. 
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Physically, the addition of heat to a supersonic flow causes it to slow down. If the thermal choking FER 
(. FERtc ) is applied, then we will have M3 = 1 (i.e. sonic combustor exit). When thermal choking occurs, it 
is not possible to increase the air mass flow through the engine. Propulsion engineers want to operate near 
thermal choking for engine efficiency reasons 5 . However, if additional heat is added, the upstream conditions 
can be altered and it is possible that this may lead to engine unstart. This is highly undesirable. For this 
reason, operating near thermal choking has been described by some propulsion engineers as “operating near 
the edge of a cliff.” In general, thermal choking will occur if FER is too high, M 0 0 is too low, altitude is too 
high too low), FTA is too high. See discussion below. 

Internal Nozzle. The exit properties M e = M e (M 3 ,A n , 7), T e = T e (M 3 . M e , 7), p e = p e (M 3 , M e , 7) of the 
scramjet’s isentropic internal nozzle are founds as follows: 


[1- 


2 ej 


M? 


T e 


Pe 


At 


t 3 


7+1 

[1 + ^ M 3 2 ] ^ 
Mf 

- 1 + 1( 7 — l)Mf - 
.1 + 1(7 — 1)M|. 


P 3 


- 1+1(7— 1 )M§ - 
.1+ 1(7-1 )M|. 


7 

7-1 


(29) 

(30) 

(31) 


where A n is the internal nozzle area ratio (see Figure 2). A n = 6.35 is used in the model. 

Thrust due to Internal Nozzle. The purpose of the expanding internal nozzle is to recover most of the 
potential energy associated with the compressed (high pressure) supersonic flow. The thrust produced by 
the scramjet’s internal nozzle is given by 47 

ThrUStinternal = r7l a (v e Vqq) + (jp e p QO )+ e (32) 

where rh a is the air mass flow through the engine, v e is the exit flow velocity, Vqo is the free-stream flow 
velocity. p e is the pressure at the engine exit plane, A\ is the engine inlet area, A e is the engine exit area, 
v e = M e sos e , t+o = MqoSOSoo, sos e = \Z^RT ei soSoo = +7 RToo , and R is the gas constant for air. Because 
we assume that the internal nozzle to be symmetric, this internal thrust is always directed along the vehicle’s 
body axis. The mass air flow into the inlet is given as follows: 


m a 


Li t a n ^ Tll ) ^ + hiCos(a) 

PooM c 


RTvc 


sin(9 s )cos(Tu ) 
sin(0 3 —a—Tu ) 


Poo d+c R T^hi COS ( 7”i i ) 


Oblique bow shock (swallowed by engine) 
Oblique bow shock - shock on lip 
Lower forebody expansion fan 


(33) 


External Nozzle. The purpose of the expanding external nozzle is recover the rest of the potential energy 
associated with the compressed supersonic flow. A nozzle that is too short would not be long enough to 
recover the stored potential energy. In such a case, the nozzle’s exit pressure would be larger than the free 
stream pressure and we say that it is under-expanded [47, Page-209]. The result is reduced thrust. A nozzle 
that is too long would result in the nozzle’s exit pressure being smaller than the free stream pressure and we 
say that it is over-expanded [47, Page-209]. The result, again, is reduced thrust. When the nozzle length is 
“properly selected,” the exit pressure is equal to the free stream pressure and maximum thrust is produced. 
Within [61, page 5], 62 the authors say that the optimum nozzle length is about 7 throat heights. This includes 
the internal as well as the external nozzle. For our vehicle, the internal nozzle has no assigned length. This 
becomes an issue when internal losses are addressed. For the Bolender, et. al. model, the external nozzle 
length is 10.15 throat heights (with throat height hi = 3.25 ft). For the new engine design presented later 
on in this paper, the external nozzle length is 7.33 throat heights (with throat height hi = 4.5 ft). The 
external nozzle contributes a force on the upper aft body. This force can be resolved into 2 components - the 
component along the fuselage water line is said to contribute to the total thrust. This component is given 
by the expression: 


ThrUStexternal 


Poo^ a 


Pe 

Poo 


In 


te) 


Pe 

Poo 


- 1 


tan(r 2 + t 1v ). 


(34) 
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Plume Assumption. The engine’s exhaust is bounded above by the aft body/nozzle and below by the 
shear layer between the gas and the free stream atmosphere. The two boundaries define the shape of the 
external nozzle. Within [60, page 1315], 7 , a critical assumption is made regarding the shape of the external 
nozzle-and-plume in order to facilitate (i.e. speed up) the calculation of the aft body pressure distribution. In 
short, the so-called “plume assumption” implies that the external nozzle-and-plume shape does not change 
with respect to the vehicle’s body axes. This implies that the plume shape is independent of the flight con- 
dition. Our (limited) studies to date show that this assumption is suitable for preliminary trade studies but 
a higher fidelity aft body pressure distribution calculation is needed to understand how properties change 
over the trimmable region. In short, our fairly limited studies suggest that the plume assumption impacts 
static properties significantly while dynamic properties are only mildly impacted. The impact of the plume 
assumption will be examined further in a subsequent publication. 

Total Thrust. The total thrust is obtained by adding the thrust due to the internal and external nozzles. 

Trimmable Region and Vehicle Properties. Within this paper (and all our work to date), trim refers 
to a non-accelerating state; i.e. no translational or rotational acceleration. Moreover, all trim analysis has 
focused on level flight. Figure 3 shows the level- flight trimmable region for the nominal vehicle being con- 
sidered'’ 19,26,29,41 (using the original nominal engine parameters). We are interested in how the static and 
dynamic properties of the vehicle vary across this region. Static properties of interest include: trim controls 
(FER and elevator), internal engine variables (e.g. temperature and pressure), thrust, thrust margin, AOA, 
L/D. Dynamic properties of interest include: vehicle instability and RHP transmission zero associated with 
FPA. Understanding how these properties vary over the trimmable region is critical for designing a robust 
nonlinear (gain-schedulted/adaptive) control system that will enable flexible operation. For example, con- 
sider a TSTO flight. The mated vehicles might fly up along q = 2000 psf to a desired altitude, then conduct 
a pull-up maneuver to reach a suitable staging altitude. 



Figure 3. Visualization of Trimmable Region: Level-Flight, Unsteady-Viscous Flow, Flexible Vehicle, 2 Controls 


III. Engine Parameter Studies 

This section examines the impact of varying the engine inlet height hi and the diffuser area ratio Ad- Three 
basic engine designs were considered: (1) current (nominal, slow or small), (2) new (intermediate speed or 
size), and (3) aggressive (fast or large). In what follows, h e denotes the internal nozzle exit height and A n 
is the internal nozzle area ratio. 
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1. Current (Nominal, Slow or Small) Engine Design. The current (nominal, slow or small) engine 
design parameters are as follows 7 : 

• hi = 3.25 h e =5 A d = 1 A n = 6.35. 

These parameters are not geometrically compatible with the vehicle shown in Figures 1 and 2; i.e. it 
would be impossible for the vehicle to have the pictorially implied flat base; i.e. internal nozzle exit 
height h e equal to inlet height h,. 

Given the above, we set out to examine engines with h e = hi. This implies that A n = -j-. 

2. New (Intermediate Speed or Size) Engine Design. The new (intermediate speed or size) engine 
design parameters were selected as follows: 

• h e = hi = 4.5 A d = 0.15 A n = A- d = 6.67. 

It should be noted that the value A d = 0.1 was used within 30,31 . This new engine design will be used 
later in the paper for analysis and control system design purposes. 

3. Aggressive (Fast or Large) Engine Design. An aggressive (fast or large) engine design was also 
considered: 

• h e = hi = 6 A d = 0.125 A n = ^ = 8. 

Constraints for Engine Parameter Trade Studies (Mach 8, 85 kft, Level Flight). The above 
engines were obtained by conducting parametric trade studies at Mach 8, 85 kft, level flight. The following 
constraints were assumed in our studies: 

• Flat base (internal nozzle exhaust height h e equal to inlet height hi); i.e. h e = h d and A n = A^ 1 ; 

• Inlet height hi was varied between ±50% of nominal 3.25 ft; 

• Engine mass rn eng i ne was varied between ±50% of nominal 10 klbs; 

• Diffuser area ratio A d was varied between 0.1 and 0.35. 

III..1. Impact of Engine Parameters on Static Properties (Mach 8, 85 kft, Level Flight) 

Figure 4 shows the impact of varying (hi, A d ) on FER, combustor temperature (assuming calorically perfect 
air), thrust, thrust margin at Mach 8, 85 kft, level flight. 

Trim FER. From Figure 4 (upper left), one observes that the: 

• trim FER decreases with decreasing A d for a fixed hp. 

• trim FER decreases with increasing h t when hi < 7. 

These suggests choosing A d small (i.e. significant diffuser compression) and h, large (i.e. large air mass flow) 
in order to achieve a small trim FER. The above, however, does not tell the full story since fuel consumption 
(trim fuel rate) - shown in Figure 4 (upper right) - increases with increasing h ,, and the thrust margin 
decreases for A d < 0.125. 

Trim Combustor Temperature. From Figure 4 (lower left), one also observes that: 

• Trim combustor temperature is a concave up function of (hi, A d ) - minimized at hi « 5.5, A d ss 0.125. 

• Trim combustor temperature exhibits a steep gradient for A d > 0.2 
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FER at Mach 8, 85kft 
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Figure 4. Trim FER., Combustor Temperature, Thrust, Thrust Margin: Dependence on hi , A a (Mach 8, 85 kft) 


Since air is assumed to be calorically perfect, it follows that high temperature effects 63 are not captured 
within the model. As such, the combustor temperatures in Figure 4 (lower left) may be excessively large. 
Future work will consider high temperature gas effects within the combustor. This is important because 
material temperature limits within the combustor are stated as 4500°i? within 64 . 

Trim Thrust Margin. From Figure 4 (lower right), we also observe that 

• Trim thrust margin is a concave down function of (hi, Ad) - maximized at hi « 6, Ad « 0.125. 

Trim Elevator and AOA. Figure 5 shows how trim elevator and AOA depend on (hi, Ad). From Figure 5, 
one observes that the: 

Elevator Deflection (deg) at Mach 8, 85kft Angle of Attack (deg) at Mach 8, 85kft 




Figure 5. Trim Elevator Deflection and Trim AOA: Dependence on ( hi , A d) - Mach 8, 85 kft, Level Flight 


• Trim elevator increases with increasing hi for a fixed Ad', 

• Trim elevator increases with decreasing Ad for a fixed hf, 
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• Trim AO A increases with increasing hi for fixed Ad- Trim AOA decreases with increasing Ad for fixed 
hi. (For hi sufficiently large, trim AOA becomes nearly independent of Ad-) 


III.. 2. Impact of Engine Parameters on Dynamic Properties (Mach 8, 85 kft, Level Flight) 


The following figure shows the impact of hi and Ad 
associated with FPA. 


Right Half Plane Pole at Mach 8, 85kft 



on the vehicle instability and RHP transmission zero 


Right Half Plane Zero at Mach 8, 85kft 



Figure 6. Right Half Plane Pole and Zero: Dependence on {hi, Ad) - Mach 8, 85 kft, Level Flight 


From Figure 6, one observes that the: 

• RHP pole increases with increasing Ad (for a fixed hi) and decreasing hi (for a fixed Ad)\ 

• RHP zero is constant with respect to Ad (for a fixed hf)\ it decreases with increasing hi (for a fixed 
Aa). 

III.. 3. Comparison of Engine Designs (Mach 8, 85 kft, Level Flight) 

In the previous sections, we considered the impact of increasing the engine height hi and diffuser area ratio 
Ad- We consider hi < 6 (bound chosen due to combustor temperature effects) and Ad > 0.125 (bound chosen 
due to thrust margin effects). Within this range, we observe the following trade-offs: 

• Increasing hi (fixed Ad) 

— PROS: Trim FER reduces, trim combustor temperature decreases (till hi ss 5.5 at Ad = 0.125), 
trim thrust margin increases, trim lift-to-drag increases (for hi > 4.0 at Ad = 0.125, not 
shown), trim drag decreases (for hi > 4.0 at Ad = 0.125, not shown), RHP pole reduces; 

— CONS: Trim fuel rate increases, trim elevator increases, trim AOA increases, RHP zero decreases, 
trim lift-to-drag decreases (for hi < 4.0 at Ad = 0.125, not shown), trim drag increases (for 
hi < 4.0 at Ad = 0.125, not shown); 

• Decreasing Ad (fixed hi) 

— PROS: Trim FER decreases, trim fuel rate decreases, trim combustor temperature decreases, trim 
thrust margin increases, RHP pole decreases (marginally); 

— CONS: Trim elevator increases, trim AOA increases (marginally), trim lift-to-drag decreases (not 
shown), trim drag increases (not shown). 

Table 4 shows a comparison of the three engine designs described above. The first is the nominal engine 
design presented in 7-9, 19 - 26 - 28 - 36 . 38 As stated earlier, this configuration is geometrically unfeasible with re- 
spect to the implied flat base vehicle diagram shown in Figures 1 and 2. As can be seen from the table, it 
is generally “slow” with a small maximum acceleration capability. The second engine design will be used 
throughout the remainder of this paper. It satisfies each of the constraints listed at the beginning of Sec- 
tion III. The third configuration is a faster configuration that also obeys the constraints. 
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Engine 

hi 

Ad 

A n 

H e 

Trim L2D 

Trim Fuel Rate 

RHP Pole 

RHP Zero 

Z/P Ratio 

Nominal 

3.25 ft 

1 

6.35 

5 ft 

2.17 

0.051 slugs/s 

3.1 

8.5 

2.7 

New 

4.5 ft 

0.15 

6.67 

4.5 ft 

4.32 

0.071 slugs/s 

2.9 

6.9 

2.4 

Fast 

6 ft 

0.125 

8 

6 ft 

4.53 

0.094 slugs/s 

2.6 

6.1 

2.3 


Engine 

Trim FER 

Trim Temp. 

Trim Thrust 

Trim Elev. 

Trim AOA 

Max Thrust 

Max Acceleration 

Nominal 

mm. 

4500°R 

1250 lbf 

9.7° 

1° 

2834 lbf 

11.1 ft/s 2 

New 


4750°R 

1380 lbf 

6.25° 

3.1° 

4647 lbf 

18.2 ft/s 2 

Fast 

0.3 

4400°R 

1445 lbf 

6.9° 

4.6° 

6582 lbf 

28.8 ft/s 2 


Table 4. Comparison of 3 Engine Designs (Mach 8, 85 kft, Level Flight) 


Table 4 shows that with respect to the nominal (slow or small) engine, the new (intermediately fast and 

sized) engine has the following associated PROS and CONS at Mach 8, 85 kft, level flight: 

• PROS: smaller trim elevator, smaller trim FER, larger maximum thrust, larger thrust margin, larger 
maximum acceleration, smaller RHP pole; 

• CONS: larger engine, larger mass, larger trim thrust, larger trim combustor temperature, larger trim 
AOA, smaller RHP zero, smaller RHP zero-pole ratio. 

Future work will examine the above tradeoffs more precisely. 

IV. Vehicle Parameter Studies 

This section examines the impact of the following vehicle parameters: 

• Engine location with respect to vehicle nose, lower forebody compression ramp inclination, center of 
gravity (CG) location with respect to vehicle nose, and vehicle mass 

across the vehicle’s level-flight trimmable region. 

IV. A. Engine Location with respect to Vehicle Nose 

In this section, we examine the impact of moving the engine rearward with respect to the vehicle nose. The 

following assumptions will be made: 

• The new engine parameters are being used; i.e. h e = hi = 4.5 ft, Ad = 0.15, A n = = 6.67. 

• The engine location is varied between 40 and 60 ft from the vehicle nose. 

• The vehicle’s center of gravity (eg) moves with the engine location (assumed fixed near engine com- 
bustor) . 

IV.A.l. Impact of Engine Location on Static Properties (Level Flight) 

Figure 7 shows how the trimmable region changes with the engine location. Specifically, it shows that the 

trimmable region shrinks as the engine is shifted rearward. 

Trim AOA. From Figure 8, we observe the following: 

• Trim AOA decreases as the engine is moved rearward for a fixed Mach and altitude. 

• For a fixed altitude, trim AOA dependence on engine location is a bit complex. When the engine is 
closer to the nose, trim AOA increases with Mach. This is because a forward situated engine results in 
a forward CG shift and hence a more stable aircraft. When the engine is moved rearward, trim AOA 
becomes insensitive to Mach variations. This dependence requires further examination. 

• For a fixed engine location and Mach, trim AOA increases with increasing altitude. 
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Figure 7. Impact of Engine Location on (Level Flight) Trimmable Region 


AOA vs. Engine Position, h=100 kft 
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Figure 8. Impact of Engine Location on Trim AOA 


smaller AOA is typically desirable when a designer wishes to increase the vehicle’s stall margin. (Generally 
eaking, an AOA that is too large results in flow separation and loss of lift.) 


Elevator vs. Engine Position, h=100 kft Elevator vs. Engine Position, M=8 



Figure 9. Impact of Engine Location on Trim Elevator 
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Trim Elevator. From Figure 9, we observe the following: 

• The trim elevator deflection increases as the engine is moved rearward for a fixed Mach and altitude. 

• For a fixed altitude, the dependence of trim elevator on engine location is nearly linear with the slope 
decreasing with increasing Mach. When the engine is closer to the nose (more stable vehicle), trim 
elevator increases with increasing Mach. When the engine is closer to the rear (more unstable vehicle), 
trim elevator increases with decreasing Mach. From this, it follows that flight at a fixed altitude and a 
low (high) Mach requires less (more) elevator for a forward situated engine (more stable vehicle), and 
more (less) elevator for a rearward situated engine (more unstable vehicle) . 

• For a fixed Mach, the dependence of trim elevator on engine location is almost linear with the slope 
increasing slightly with increasing altitude. When the engine location is also fixed, the trim elevator 
increases with increasing altitude. 

• Trim elevator deflection increases with increasing altitude. 

Trim FER. Figure 10 illustrates how trim FER depends on engine location. From the figure, the following 



is observed: 

• Trim FER is a concave up function with respect to engine location for a fixed altitude and Mach - with 
trim FER being minimized near 45 ft for most flight conditions. 

• For a fixed engine location and altitude (or Mach), trim FER increases with increasing Mach (or 
altitude). 

IV. A. 2. Impact of Engine Location on Dynamic Properties (Level Flight) 


RHP Pole vs. Engine Position, h=100 kft RHP Pole vs. Engine Position, M= 8 



Figure 11. Impact of Engine Location on Right Half Plane Pole 


RHP Pole. Figure 11 shows that the 
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• instability increases (roughly linearly) as the engine is moved rearward; 

• instability increases with increasing Mach and decreasing altitude. 

Moving the engine rearward, moves the center of gravity (eg) rearward with respect to the aerodynamic 
center (ac) - thus making the airplane more unstable. 

Motivation for Increased Instability. From the above, a designer may wish to increase the vehicle insta- 
bility in order to make the vehicle more maneuverable in terms of following aggressive flight path angle 
or vertical acceleration commands. It may also be desirable in order to facilitate the attenuation of high 
frequency wind disturbances. (The link between instability and maneuverability was understood by the 
Wright Brothers early on in their work [65, page 39]. This could be important for a missile going after agile 
targets. Such might be the case for military applications. In such a case, one should note that a larger 
instability requires a larger minimum control system bandwidth for vehicle stabilization. 41 This, however, 
may conflict with higher frequency non-minimum phase, structural, aero-elastic, and actuator dynamics. In 
the same spirit, a larger bandwidth at the elevator would typically require a faster control surface actuator. 
Such considerations must be rigorously addressed at some point in the design process - the sooner, the better. 


RHP Zero vs. Engine Position, h=100 kft RHP Zero vs. Engine Position, M= 8 



Figure 12. Impact of Engine Location on Right Half Plane Zero 


RHP Zero. From Figure 12, one observes that the: 

• RHP zero varies little with engine position for a fixed altitude and Mach 

• RHP zero increases with increasing Mach and decreasing altitude 

Increasing RHP Zero: Moving Engine Rearward. From the above, it follows that one might move the engine 
rearward (making the vehicle more unstable) in order to maximize the right half plane zero. By so doing, 
a vehicle designer can (in principle) increase the maximum achievable flight path angle bandwidth. 41 One 
must, of course, note that flexible modes, the associated uncertainty, and the control system simplicity can 
also limit the achievable bandwidth. This will be the case when the flexible modes lie within a decade of the 
right half plane zero. Additional pros associated with moving the engine rearward include: less fuel usage 
- minimized near 55 ft (not shown). Associated cons include the following: trim L/D drops (monotonically 
for 40 to 60 ft interval), trim FER increases (FER/thrust/acceleration margin decreases). 

IV.B. Lower Forebody Inclination 

In this section, we examine the impact of varying the lower forebody inclination angle. The following is 
assumed: 

• New engine parameters; i.e. h e = hi = 4.5 ft, Ad = 0.15, A n = = 6.67. 

• Lower forebody inclination varied from 4.2° to 8.2° (nominal value = 6.2°) 

• All lengths (forebody, aftbody, engine length), upper forebody angle kept constant 
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AOA (deg) 


• Tail angle and total vehicle height change as a result 

• CG assumed to be fixed 

• Heating effects due to slender nose are not considered 

IV.B.l. Impact of Lower Forebody Inclination on Static Properties (Level Flight) 

Trimmable Region. Figure 13 shows how the trimmable region changes with lower forebody inclination angle. 
From Figure 13, one observe that the: 



Mach 

Figure 13. Impact of Lower Forebody Inclination on (Level Flight) Trimmable Region 


• Trimmable region shrinks with increasing lower forebody inclination 

• Pinch point moves toward a higher Mach and lower altitude with increasing lower forebody inclination 

Trim AOA. Figure 14 shows how AOA varies with lower forebody inclination angle. From Figure 14, one 
observes that the: 



Figure 14. Impact of Lower Forebody Inclination on Trim AOA 


• Trim AOA decreases linearly with increasing lower forebody inclination for a fixed altitude 
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• Trim AOA decreases with Mach at lower Mach numbers for a fixed altitude; 


• Trim AOA increases with increasing altitude for a fixed Mach 

Trim Elevator. Figure 15 shows how elevator varies with lower forebody inclination angle. From Figure 15, 
one observes that the: 



Figure 15. Impact of Lower Forebody Inclination on Trim Elevator 


• Trim elevator deflection increases linearly with increasing forebody inclination 

• Trim elevator deflection increases with increasing Mach, increasing altitude 
Trim FER. Figure 16 shows how FER varies with lower forebody inclination angle. 



Lower forebody inclination (deg) Lower forebody inclination (deg) 


Figure 16. Impact of Lower Forebody Inclination on Trim FER 


From Figure 16, one observes that the: 

• Trim FER increases almost linearly with increasing lower forebody inclination 

IV. B. 2. Impact of Lower Forebody Inclination on Dynamic Properties (Level Flight) 

RHP Pole. Figure 17 shows how the RHP pole varies with lower forebody inclination. From Figure 17, one 
observes that the: 

• RHP pole increasing with increasing lower forebody inclination 

• RHP pole increases with Mach and decreasing altitude 

RHP Zero. Figure 18 shows how the RHP pole varies with lower forebody inclination. From Figure 18, one 
observes that the: 

• RHP zero increases with increasing lower forebody inclination 

• RHP zero increases with increasing Mach and decreasing altitude 
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RHP Pole vs. Lower forebody inclination, h=1 00 kft RHP Pole vs. Lower forebody inclination, M= 8 



Figure 17. Impact of Lower Forebody Inclination on Right Half Plane Pole 


RHP Zero vs. Lower forebody inclination, h=100 kft RHP Zero vs. Lower forebody inclination, M= 8 



Figure 18. Impact of Lower Forebody Inclination on Right Half Plane Zero 


IV. C. Center of Gravity 

This section examines the impact of varying the center of gravity (eg) location. The following assumptions 
are made: 

• New engine parameters; i.e. h e = hi = 4.5 ft, Ad = 0.15, A n = = 6.67. 

• CG varied from 45ft to 65ft artificially (no internal changes are made to achieve this shift) 

• For dynamic properties, we only consider CG locations for which vehicle is open loop unstable 

IV.C.l. Impact of Center of Gravity on Static Properties (Level Flight) 

Trimmable Region. Figure 19 shows how the trimmable region changes with a shifting vehicle CG. From 
Figure 19, one observes that the 

• trimmable region decreases slightly as the vehicle CG is moved rearward. 

Trim AOA. Figure 20 shows how trim AOA depends on the vehicle’s CG location. 

Figure 20 shows that the 

• Trim AOA decreases as the vehicle CG moves rearward. 

Trim Elevator. Figure 21 shows how trim elevator depends on eg location. 

• Trim elevator deflection increases as the vehicle CG moves rearward. 

Trim FER. Figure 22 shows how trim FER depends on vehicle CG location. Figure 22 shows that the 

• Trim FER increases as the vehicle CG moves rearward. 
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Elevator deflection (deg) AOA (deg) 


Envelope Variations with CG Location 



AOA vs. CG Location, h=100 kft AOA vs. CG Location, M=8 



Figure 20. Impact of Vehicle CG Location on Trim AOA 


Elevator vs. CG Location, h=100 kft Elevator vs. CG Location, M=8 



IV. C. 2. Impact of Center of Gravity on Dynamic Properties (Level Flight) 

RHP Pole. Figure 23 shows how the vehicle instability depends on the CG location. Figure 23 shows that 
the vehicle 
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FER vs. CG Location, h=1 00 kft FER vs. CG Location, M=8 



Figure 22. Impact of Vehicle CG Location on Trim FER, 


RFIP Pole vs. CG Location, h=1 00 kft RFIP Pole vs. CG Location, M= 8 



Figure 23. Impact of Vehicle CG Location on Right Half Plane Pole 


• RHP pole increases linearly as the vehicle CG moves rearward. 


RHP Zero. Figure 24 shows how the RHP zero depends on the CG location. Figure 24 shows that the 
vehicle 


RFIP Zero vs. CG Location, h=100 kft 



RFIP Zero vs. CG Location, M= 8 
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Figure 24. Impact of Vehicle CG Location on Right Half Plane Zero 


• RHP zero decreases linearly as the vehicle CG moves rearward. 

IV. D. Vehicle Mass 

This section examines the impact of varying the vehicle’s (total) mass. The following assumptions are made: 

• New engine parameters; i.e. h e = hi = 4.5 ft, Ad = 0.15, A n = = 6.67. 
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AOA (deg) 


• Mass of vehicle modified without changing any material or subsystem properties. 
IV.D.l. Impact of Vehicle Mass on Static Properties 

Trimmable Region. Figure 25 shows how the trimmable region depends on vehicle mass. 


Envelope Variations with Mass Ratio 



Figure 25. Impact of Vehicle Mass on (Level Flight) Trimmable Region 


From Figure 25, one observes that the: 

• Trimmable region shrinks as the vehicle mass is increased. 

Trim AOA. Figure 26 shows how trim AOA depends on vehicle mass. The figure shows that the: 


AOA vs. Mass Ratio, h=100 kft A0A vs - Mass Rati °’ M = 8 



Figure 26. Impact of Vehicle Mass on Trim AOA 


• Trim AOA increases as the vehicle mass is increased. 

Trim Elevator. Figure 27 shows how trim elevator depends on vehicle mass. The figure shows that the: 

• Trim elevator deflection increases as the vehicle mass is increased. 

Trim FER. Figure 28 shows how trim FER depends on vehicle mass. The figure shows that the: 

• Trim FER increases as the vehicle mass is increased. 
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Elevator vs. Mass Ratio, h=100 kft Elevator vs. Mass Ratio, M=8 



Figure 27. Impact of Vehicle Mass on Trim Elevator 



Figure 28. Trim FER vs Vehicle Mass 


IV. D. 2. Impact of Vehicle Mass on Dynamic Properties (Level Flight) 

RHP Pole. Figure 29 shows how the vehicle instability depends on vehicle mass. The figure shows that the: 


RHP Pole vs. Mass Ratio, h=100 kft RHP Pole vs. Mass Ratio, M= 8 



Figure 29. Impact of Vehicle Mass on Right Half Plane Pole 


• RHP pole decreases as the vehicle mass is increased. 


RHP Zero. Figure 30 shows how the RHP zero depends on vehicle mass. The figure shows that the: 


• RHP zero decreases as the vehicle mass is increased. 
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RHP Zero vs. Mass Ratio, h=1 00 kft RH p Zero vs. Mass Ratio, M= 8 



Figure 30. Impact of Vehicle Mass on Right Half Plane Zero 


IV. E. Summary of Level Flight Trim Trade Studies 

The following tables summarize each of the conducted level- flight trim trade studies. 

Rearward Engine Shift. Table 5 summarizes trends for a rearward shifting engine. The engine is shifted 
rearward (with the CG) with the vehicle height kept constant; i.e. the lower forebody inclination angle is 
decreasing, thus making the vehicle sharper. As the engine is shifted rearward, we observe specific monotonic 
trends and tradeoffs that result in the following PROS and CONS: 

• PROS: trim AOA decreases, trim lift remains nearly constant; 

• CONS: trim elevator increases (CG rearward), trim drag increases, trim L/D decreases, RHP pole 
increases. 

We also observe the following more complex (non-monotonic) behavior: 

• trim fuel rate decreases (min near 55 ft) and then increases; 

• trim FER decreases (min near 45 ft - thrust/acceleration margin increases) and then increases (thrust/acceleration 
margin decreases). 


Property 

Pro 

Con 

Trim Lift 

Almost Constant 

Almost Constant 

Trim Drag 


Increases monotonically 

Trim L/D 
Trim AOA 

Decreases monotonically 

Decreases monotonically 

Trim Elevator 


Increases monotonically 

Trim FER 

Decreases till 45 ft 

Increases after 45 ft 

Trim Fuel Rate 

Decreases till 55 ft 

Increases after 55 ft 

RHP Pole 


Increases monotonically 

RHP Zero 

Increases after 45 ft 

Decreases till 45 ft 

RHP Z/P Ratio 


Decreases monotonically 


Table 5. Trends for Rearward Engine Shift 


Lower Forebody Angle Increase. Table 6 summarizes trends for an increasing lower forebody angle. Here, 
the horizontal engine location is fixed and the engine is moved downward - thus increasing the height of the 
vehicle. As the lower forebody angle is increased, we observe specific monotonic trends and tradeoffs that 
result in the following PROS and CONS: 
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• PROS: trim AOA decreases, RHP zero increases, RHP zero increases, RHP zero-pole ratio almost 
constant; 

• CONS: trim elevator increases, trim drag increases, trim L/D decreases, trim FER increases, trim fuel 
rate increases, RHP pole increases. 

Non-monotonic patterns are not observed in this case. 


Property 

Pro 

Con 

Trim Lift 
Trim Drag 

Almost Constant 

Almost Constant 

Trim L/D 


Decreases monotonically 

Trim AOA 

Decreases monotonically 


Trim Elevator 


Increases monotonically 

Trim FER 


Increases monotonically 

Trim Fuel Rate 

Decreases monotonically 


RHP Pole 


Increases monotonically 

RHP Zero 

Increases monotonically 


RHP Z/P Ratio 

Almost constant 

Almost constant 


Table 6. Trends for Increasing Lower Forebody Inclination 


Rearward CG Shift. Table 7 summarizes trends for a rearward shifting CG. As the CG is shifted rearward, 
we observe specific monotonic trends and tradeoffs that result in the following PROS and CONS: 

• PROS: trim AOA decreases, trim fuel rate almost constant; 

• CONS: trim elevator increases, trim drag increases, trim L/D decreases, trim FER increases, RHP pole 
increases, RHP zero decreases, RHP zero-pole ration decreases. 

Non- monotonic patterns are not observed in this case. 


Property 

Pro 

Con 

Trim Lift 

Almost Constant 

Almost Constant 

Trim Drag 


Increases monotonically 

Trim L/D 
Trim AOA 

Decreases monotonically 

Decreases monotonically 

Trim Elevator 


Increases monotonically 

Trim FER 


Increases monotonically 

Trim Fuel Rate 


Almost constant 

RHP Pole 


Increases monotonically 

RHP Zero 


Decreases monotonically 

RHP Z/P Ratio 


Decreases monotonically 


Table 7. Trends for Rearward CG Shift 


Vehicle Mass Increase. Table 8 summarizes trends for increasing mass. As the vehicle mass is increased, we 
observe specific monotonic trends and tradeoffs that result in the following PROS and CONS: 

• PROS: RHP pole decreases, RHP zero-pole ratio almost constant; 

• CONS: trim elevator and trim AOA increase, trim FER and fuel rate increase, RHP zero decreases. 
We also observe the following more complex (non-monotonic) behavior: 
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• trim L/D is a concave down function of mass for a fixed Mach and altitude; mass at which peak occurs 
increases with Mach (altitude fixed); decreases with altitude (Mach fixed). 


Property 

Pro 

Con 

Trim Lift 

Increases monotonically 


Trim Drag 


Increases monotonically 

Trim L/D (peak) 

Increases with Mach (altitude fixed) 

Decreases with altitude (Mach fixed) 

Trim AOA 


Increases monotonically 

Trim Elevator 


Increases monotonically 

Trim FER 


Increases monotonically 

Trim Fuel Rate 


Increases monotonically 

RHP Pole 

Decreases monotonically 


RHP Zero 


Decreases monotonically 

RHP Z/P Ratio 

Almost constant 

Almost constant 


Table 8. Trends for Increasing Mass 


V. Vehicle Optimization 

Various schemes have been considered for the optimization of space vehicles. Within 66 , a conceptual 
design process - that takes factors like cost analysis into consideration - is described. The non-hierarchical 
nature of the design process is illustrated and several optimization methods - parameter based, gradient 
based, and stochastic methods - are examined. Within 67 , a probabilistic approach to vehicle design is 
taken in order to account for uncertainties. For air-breathing hypersonic aircraft, the coupling between the 
airframe and the engine introduces additional constraints into the design process. The importance of a 
multidisciplinary design optimization approach is illustrated within 68 . A collaborative approach 69 for launch 
vehicle design is considered within 70 . The design is decomposed into several subsystem design problems. 
A system level optimizer coordinates the integration of subsystem designs while taking into account inter- 
subsystem coupling and constraints. Such a subsystem-based algorithm is well suited to exploit parallel 
computing architectures. 

In this section, we address optimizing the generic carrot shaped hypersonic vehicle under consideration. 
The flat base assumption is made; i.e. h e = hi , A n = -j-. The goal is twofold: (1) understand what config- 
urations result when very simple optimizations are considered, (2) understand how the gap matric can be 
used to design a control- friendly vehicle. 


Admissible Parameter Values. Optimizations are conducted over the following control surface and engine 
parameter values: 


• Elevator area S e i ev : 

• Horizontal elevator position x e i e 

• Inlet height hf. 

• Lower forebody inclination til- 


[8.5, 25.5] ft 2 (nominal: 17 ft 2 ) 

[70, 90] ft from nose (nominal: 85 ft) 
[2.25, 6.75] ft (nominal: 4.5 ft) 

[4.2°, 8.2°] (nominal: 6.2°) 


Overview of Trim Optimization Algorithm. A nonlinear simplex optimization method was used. No 
derivatives are calculated because trimming the vehicle results in each iteration being expensive. 71, 72 An 
overview of the vehicle optimization algorithm is as follows: 

1. Choose cost function (e.g. trim AOA, trim fuel, trim elevator, trim FER, trim L/D); 

2. Select values for parameters being optimized over; e.g. S e i ev , x e i ev , hi, til;, 

3. Trim the vehicle at selected flight condition; i.e. Mach 8, 85 kft, level flight; 
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4. Evaluate cost function using obtained trim values; 

5. Repeat 2-4 until optimal configuration for chosen cost function is obtained. 

ASU High performance Computing Cluster. All optimizations presented in this paper were performed 
using the Arizona State University (ASU) High Performance Computing (HPC) cluster. This cluster is part 
of ASU’s Ira A. Fulton School of Engineering High Performance Computing Initiative. The HPC is composed 
of 220 dual quadcore Intel Xeon EM64T (Extended Memory 64 Technology) nodes - each with 16 gigabytes 
of RAM. The cluster also has a partition for running large numbers of serial jobs comprised of 185 nodes 
with dual Xeon MP 64bit processors. 

Static Optimization Using ASU HPC: Computational Issues. While individual optimizations pro- 
ceed in a serial fashion, the search space is partitioned into subspaces. Due to the nonlinear nature of 
the problem, a multistart algorithm was used (with at most 4 4 = 256 initial conditions per optimization 
problem); i.e. several initial starting points are considered in each subspace. Optimization from an initial 
guess for a single operating point using the nonlinear (3DOF + flexibility) model takes approximately 15 
minutes (on a 2.5GHz CPU with 1GB of memory) to converge to a local minimum (for the static objectives). 
The results from all the initial conditions are examined and the optimum is chosen. On average, 256 initial 
guesses are chosen for each static objective minimization. 

The q = 2000 psf fuel consumption objective optimization takes approximately an hour to converge to a 
local minimum (on a 2.5GHz CPU with 1GB of memory) from an initial guess (in each iteration the plant 
must be trimmed at 7 points along the q 2 ooo trajectory). Similar to the static optimization, a multistart 
algorithm was used (with at most 4 4 = 256 initial conditions) . 

Trim Optimization Results (Mach 8, 85kft, Level Flight). Table 9 shows the results obtained by 
optimizing the following trim vehicle variables at Mach 8, 85kt (level flight): trim AOA, trim fuel, trim 
elevator deflection, trim FER, and trim L/D. Table 10 shows the parameters that result from the respective 
optimizations. Table 11 contains corresponding dynamical characteristics. 


Objective 

Trim 

L/D 

Trim 

FER 

Trim Fuel 
Rate(^f) 

Trim 

Elevator (deg) 

Trim 

AOA (deg) 

(1) Maximize Trim L/D 

6.63 

0.29 

0.116 

2.19 

6.39 

(2) Minimize Trim FER 

5.95 

0.27 

0.102 

2.65 

6.64 

(3) Minimize Trim Elevator 

5.95 

0.27 

0.102 

1.11 

6.64 

(4) Minimize Fuel Consumed (const, q) 

5.89 

0.27 

0.102 

3.11 

6.64 

(5) Minimize Trim Fuel Rate 

4.38 

0.47 

0.0361 

3.93 

2.19 

(6) Nominal (with New Engine) 

4.30 

0.35 

0.0708 

6.44 

2.88 

(7) Minimize Trim AOA 

3.28 

0.49 

0.0616 

10.14 

0.00271 


Table 9. Optimized Static Properties - Resulting Trim Values (Mach 8, 85 kft, Level Flight) 


Objective 

Selev (ft 2 ) 

%elev 

hi (ft) 

til (deg) 

(1) Maximize Trim L/D 

25.5 

90 

6.75 

5.24 

(2) Minimize Trim FER 

19.55 

90 

6.75 

422 

(3) Minimize Trim Elevator 

25.5 

90 

6.75 

422 

(4) Minimize Fuel Consumed (const, q) 

18.19 

90 

6.75 

422 

(5) Minimize Trim Fuel Rate 

16.15 

81.76 

2.25 

422 

(6) Nominal (with New Engine) 

17 

85 

4.5 

6.2 

(7) Minimize Trim AOA 

19.89 

73 

3.15 

O 


Table 10. Optimized Static Properties - Resulting Vehicle Configurations/Parameters (Mach 8, 85 kft, Level Flight) 
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Objective 

RHP Pole 

RHP Zero 

RHP Zero-Pole ratio 

(1) Maximize Trim L/D 

(2) Minimize Trim FER 

(3) Minimize Trim Elevator 

(4) Minimize Fuel Consumed (const, q) 

(5) Minimize Trim Fuel Rate 

1.60 (-45.18%) 
1.95 (-33.43%) 
1.48 (-49.56%) 
2.05 (-30.07%) 
3.14 (+7.23%) 

5.83 (-16.02%) 
5.62 (-19.06%) 
5.80 (-16.46%) 
5.61 (-19.15%) 
7.32 (+5.44%) 

3.6325 (+53.18%) 
2.8833 (+21.59%) 
3.928 (+65.64%) 
2.7416 (15.61%) 
2.3319 (-1.66%) 

(6) Nominal (with New Engine) 

2.93 (0.00%) 

6.94 (0.00%) 

2.37 (0.00%) 

(7) Minimize Trim AOA 

3.55 (+21.23%) 

6.34 (-8.59%) 

1.788 (-24.60%) 


Table 11. Optimized Static Properties - Resulting Dynamic Characteristics (Mach 8, 85 kft, Level Flight) 


Within Table 10, over bars imply that the maximum admissible parameter was achieved; underlines imply 
that the minimum admissible parameter bound was achieved. From Tables 9-10, one observes the following: 


1. Maximizing Trim L/D. Maximizing trim L/D results in a small FER - marginally above that obtained 
when trim FER is directly minimized. As when trim FER is minimized, maximizing L/D also results in 
a large trim AOA. In short, maximizing trim L/D yields results comparable to minimizing trim FER. 
Maximizing trim L/D, however, results in a larger elevator surface area and lower forebody inclination 
vis-a-vis when trim FER is minimized. For this case, S e i ev , x e i ev , and hi, were maximized to yield a 
maximally large elevator, a maximally effective elevator, a maximally large engine. 


2. Minimizing Trim FER. Minimizing trim FER produces results that are comparable to those obtained 
when L/D is maximized - with L/D reduced 10.3% from the maximum L/D, FER reduced by 12.1%, 
and AOA increased by only 3.9%. In contrast to maximizing L/D, minimizing FER results in a much 
smaller elevator surface area (23.3%) and lower forebody angle (19.9%). For this case, x e i ev , and hi 
were maximized while t\l was minimized to yield a maximally effective elevator, a maximally large 
engine, and a very aerodynamic lower forebody. 

3. Minimizing Trim Elevator. Minimizing trim elevator produces results (L/D, Fuel, FER, AOA) that 
are comparable to minimizing trim FER with less than half the elevator. For this case, S e i ev , x e i ev , and 
hi were maximized while til was minimized to yield a maximally large elevator and effective elevator, 
a maximally large engine, and a very aerodynamic lower forebody. 

4. Minimizing Fuel Consumed (const, q). Minimizing the fuel consumed along a constant dynamic 
pressure trajectory produces results that are comparable to minimizing trim FER with x e i ev and hi 
maximized and til, minimized to yield a maximally effective elevator, a maximally large engine, and a 
very aerodynamic lower forebody. This optimization is more involved than the others being considered. 
Hence some additional explanation is required. 

The objective here is to minimize the total fuel consumed while flying along the q = 2000 psf altitude- 
Mach profile from Mach 5.52 at 70kft to Mach 11.08 at a lOOkft. To do so, we approximate the total 
fuel consumed. For simplicity, vehicle mass changes are ignored. The method used to approximate fuel 
consumption is described below. 


(a) Trim at Selected Points. Select points along the stated q = 2000 psf profile spaced 5 kft. This 
yields the 

altitudes: [70, 75, 80, 85, 90, 95, 100] kft and the corresponding 

Machs: [5.52, 6.21, 6.98, 7.85, 8.81, 9.88, 11.08]. 

This divides the profile into 6 legs. Trim the vehicle at each point. 


(b) Approximate Maximum Accelerations. Approximate the maximum horizontal acceleration possi- 
ble at the left end point of each leg. Do so using the FER margin at each left end point. The FER 
margin is used to determine the maximum FER. This is used to compute the associated thrust 
Th(FER max )i. The trim drag Di and trim AOA cq are then used to approximate the maximum 


acceleration at the left end point: at 
vehicle. 


Th(FER rnax )i cos oti-Dj 


where m is the total mass of the 
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(c) Determine Average Accelerations for Each Leg. Assume that the average acceleration in leg i 
is the mean of the two end-point accelerations; i.e. a, = a *~ l ~°»+ 1 where is the approximate 
maximum acceleration at the left end point of the z th leg and a *+ 1 is that at the right end point. 

(d) Approximate Time to Fly Each Leg. Approximate the time to fly leg i as T) « t, * + ir Ut . 

(e) Approximate the average fuel rate for leg i (P,) as the mean of the fuel rates at the end points 
(where we have maximum acceleration); i.e. P, = F,+F,+1 

(f) Approximate Fuel Consumed During Each Leg. Approximate the fuel consumed during leg i as 

Fi ~ Fi x Ti 

(g) Approximate Total Fuel Consumed Over Flight Profile. Approximate the total fuel consumed 
(objective to be minimized) as 


6 

Total fuel consumed ss F t (35) 

i - 1 


(h) Inadmissible Designs. We also want to maintain a reasonable FPA. The average FPA for the i th 
leg is approximated as 


hi - (-1 hi 

Vi x Ti 


(36) 


where (u,, hi) are the velocity and altitude at the left end point of the i th leg and hi+\ is the 
altitude at the right end point. Designs for which the average FPA in any leg exceeded 3° were 
deemed inadmissible and discarded. 


5. Minimizing Trim Fuel Rate. Minimizing trim fuel rate produces results that lie between those obtained 
when trim AOA is minimized and elevator is minimized (or trim FER is minimized or trim L/D is 
maximized). Surprisingly, for this case both hi and t\l are minimized to yield a maximally small 
engine with a very aerodynamic lower forebody. This optimization results in low trim air mass flow 
to the engine (and low trim AOA) so that even with the resulting larger trim FER and increased trim 
thrust, the trim fuel rate is small. 

6. Nominal Vehicle. Observe that the nominal values S e i ev = 17, hi = 4.5, t\l = 6.2 are midpoints 
while the nominal value x e i ev = 85 is near the right end point. Note that this case can be viewed as 
a tradeoff between minimizing trim fuel rate and trim AOA when looking at trim L/D and elevator, 
elevator area, and lower forebody angle. Otherwise, trim FER, fuel rate, AOA, elevator rearward 
distance, and engine inlet height are larger. 

7. Minimizing Trim AOA. Minimizing trim AOA uses excessive trim elevator for lift. This results in 
excessive drag and results in a small trim L/D. It also produces a large trim FER and small trim fuel 
rate. The resulting low trim fuel rate is due to the reduced trim air mass flow to the engine (due to 
the smaller trim AOA). From a vehicle configuration perspective, this case results in the largest lower 
forebody inclination (least aerodynamic shape) - achieving the maximum admissible value for t\l- 

Motivation for Gap Metric Based Optimization. The gap metric represents a system-theoretic measure 
that quantifies the “distance” between two dynamical systems and whether or not a common controller can 
be deployed for the systems under consideration 48 . Within 49 , the gap between two LTI dynamical systems 
(Pi, P 2 ) is defined as follows: 

ff (Pi, P 2 ) =‘ max{ Q mf^ II [ ]-[£“] Q IU Q |nf^ II [ n\ ] ~ [ n\ ] Q lloo } (37) 

where Pi = NiDf 1 , P 2 = 7V 2 Df 1 , and (Ni,Di) denotes a normalized right coprime factorization for L\ 
(i = 1,2) in the sense of 73 . The gap metric (and the v gap 74 ) has often been considered from a robustness 
perspective in the stabilization of feedback systems' 5 . Within 76 , the authors relate the gap metric with tra- 
ditional stability margins. The gap metric has also been considered for the design of controllers for space 
vehicles 77, ' 8 . Given the importance of this metric from a controller design perspective, it may be very useful 
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to consider the gap metric during the vehicle design phase. This motivates the vehicle optimization which 
we consider below. 

Vehicle Optimization Using Gap Metric: Minimizing Plant Variability. We now consider vehicle 
optimization using the gap metric. Details of the optimization are now given. 

• The optimization is conducted with respect to the four parameters: S e i ev , x e i ev , hi, til- 



Mach 

Figure 31. Directed Set of Level Flight Operating Points for Cumulative Gap Calculation 


• Figure 31 shows the set of operating points used - a directed set - over which a cumulative gap is 
defined. The figure shows 17 operating points. (The operating points were selected in a way so as 
to span the desired trimmable region. How to do this systematically will be pursued in future work.) 
Associated with these 17 points are 16 gaps: gapi = gap(Pi, Pi+i), i = 1,2,. ..,16. The cumulative 
gap associated with these 17 points (along the directed set) is defined as follows: 

16 

cumulative gap = £ gapi (38) 

*=i 

A vehicle configuration is selected during optimization. The vehicle is trimmed at each point along the 
directed set. The gap between each pair of neighbors (along the directed set) is calculated and added 
to the running total for the vehicle configuration being considered. 

• If a particular vehicle configuration is not able to trim at some point along the directed set, a penalty 
is added to the cumulative gap in lieu of the gap between this frozen point and its neighbors (along 
the directed set). 

Gap Optimization Using ASU HPC: Computational Issues. Obtaining a gap-optimized vehicle from 
an initial guess takes approximately 4 hours on an ASU HPC 2.5GHz CPU with 1GB of memory. (During 
each iteration, the plant must be trimmed and linearized at each of the 17 selected operating points in 
altitude-Mach space.) As for the static optimization, a multistart algorithm was used (with at most 3 4 = 81 
initial conditions). 

Vehicle Optimization Results: Gap-Optimized Versus Nominal. The parameters for the nominal 
vehicle (with new engine design) are as follows: 

• Nominal Vehicle Parameters (with New Engine) 

Seiev = PI ft 2 X e i ev = 85 ft hi = 4.5 ft Til = 6.2° 

The cumulative gap for the nominal vehicle (with new engine design) is 4.4994 (mean gap = 0.28121). 
The final parameters for the gap-optimized vehicle are as follows: 

• Gap-Optimized Vehicle Parameters 

S elev = 25.5 ft 2 x elev = 90 ft hi = 6.12 ft r 1L = 4.2° 

The cumulative gap for the optimized vehicle is 3.3989 (mean gap = 0.21243). 
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Note that with respect to the nominal vehicle, the gap-optimized vehicle has a 

• larger elevator which is slightly shifted rearward (more effective elevator), larger engine inlet height 
(faster vehicle) , smaller lower forebody inclination angle (more aerodynamic) . 

Modulo the simplicity of the nonlinear (3DOF + flexibility) model being used in this study and its associ- 
ated limitations, the results given below suggest that these trends may significantly simplify control system 
design. More comprehensive studies are (of course) required in order to make definitive claims. Such a study 
could, for example, involve the 6DOF models being developed within 79 . 


RHP Pole Variation Over Level Flight Trimmable Region: Gap-Optimized Versus Nominal. 

Figure 32 (left) shows how the RHP pole varies across the trimmable region when the cumulative gap is 
minimized across the trimmable region (as discussed above). The variation for the nominal system (using 
the new engine design) is also given in Figure 32 (right). 


Gap Optimized Plant - RHP Pole 


Nominal Plant - RHP Pole 




Figure 32. RHP Pole Over Trimmable Region - Nominal and Gap Optimized Vehicle 


Figure 32 (left) shows that the 

• RHP pole variation has decreased significantly when the vehicle is gap optimized; we will see that this 
implies that less adaptation (i.e. scheduling) will be required in the final control system design with 
respect to that required for the system associated with Figure 32 (right); 

• Trimmable region has been expanded at low Mach and altitudes 


VI. Control System Design Studies 

A classic decentralized inner-outer loop control system architecture was used to illustrate control design 
issues. Such an architecture was examined within [41]. It can be visualized as shown in Figure 33. 



Figure 33. Inner Outer Loop Control System 
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Outputs: y p = [ v 7 ] T 


Controls: u = [ FER 5 e ] T 


State Feedback: x r = [ 6 d] T 


Overview of Nominal Control System Architecture. Within Figure 33, we have an inner loop con- 


troller: 



0 


0 


— 9i z i ff 

s 2 +bis+b 0 

-9i% 

s 2 +bis+b 0 


s z -\-a\s-\-a 0 

s z -\-a\s-\-a 0 
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(39) 


and an outer loop controller: 


K 0 (s) 


g b (s+z b ) 
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10 

s+10 


1 3 


0 


g^{s-\-z^ 


20 


s+20 


(40) 


The inner-outer loop structure is now described. 

1. Velocity Control Via Single Loop PI. A PI controller (with roll off) is used to control velocity. Two 
parameters are associated with velocity control: ( gt , Zb ). A simple anti- windup method (not discussed 
in the paper) is used to address the state-dependent thermal choking nonlinearity discussed earlier. 

2. FPA Control Via Inner-Outer Loop. A PD (proportional plus derivative with notch) inner-loop con- 
troller is used on pitch to stabilize the vehicle’s pitch dynamics and make the modified dynamics look 
friendly to the outer-loop FPA PI (proportional plus integral) controller. The inner loop’s notch pro- 
vides lead to assist with the first flexible mode’s complex zero-pole pair (near 21 rad/sec). As might 
be expected with this inner-outer loop structure, the inner loop is generally much faster than the outer 
loop. Its bandwidth must be large enough to stabilize the vehicle and follow outer-loop commands, but 
not too high so that the flexible dynamics are overly excited. (It should be noted that any practical 
“real-world” implementation of the PD controller should include additional roll off in order to attenuate 
high frequency noise n, resulting from 9 and theta measurements.) 

A PI (proportional plus integral with roll off) outer-loop controller is used for FPA. The outer loop is 
generally much slower than the inner loop. The maximum achievable outer loop bandwidth is limited 
by the vehicle’s flexible dynamics as well as the right half plane zero associated with the elevator-FPA 
map. 

Four parameters are associated with FPA control: (gi , Zi) for the inner-loop PD controller, g 1 , z 7 ) for 
the outer-loop PI controller. 

Command pre-filters are also included: \\h, = •— for velocity reference commands, Wfpa = V, for FPA 

reference commands. (These reference command pre-filters are not shown within Figure 33.) 

Additional information about controller parameter selection can be found within 41 . Nominal control 
system parameter values are given in Table 12. 


9i 

Zi 

9j 

z i 

9b 

Zb 

o 0 

ai 

bo 

bi 

1 

4.5 

21.23 

0.404 

2.2437 

0.0278 

30 

900 

15 

250 


Table 12. Control System Parameters for Nominal Hypersonic Vehicle (Mach 8, 85kft, Level Flight) 


VI. A. Impact of Engine Location on Control System Design 

This section shows how a nominal control system design is impacted by the selection of the engine location 
with respect to the nose. 

Impact on Velocity Loop’s PI. As the engine is moved rearward (and the vehicle becomes more unstable), 
the (see Figure 34 - right) 

• velocity loop’s PI controller gain increases (smaller slope than FPA PI gain) 
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Velocity PI 
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Figure 34. Control System Parameter Dependencies on Engine Location 


• velocity loop’s PI controller zero remains constant 

Impact on FPA Inner-Loop PD. As the engine is moved rearward (and the vehicle becomes more 
unstable), the 

• inner-loop PD gain remains constant 

• inner-loop PD zero increases linearly - due to linearly increasing RHP Pole (see Figure 11) 

Impact on FPA Outer-Loop PI. As the engine is moved rearward (and the vehicle becomes more 
unstable), the 

• outer-loop PI controller gain increases; as inner-loop becomes more aggressive (due to increasing PD 
zero), the outer- loop must compensate - hence similar slope on outer-loop gain 

• outer-loop PI controller zero increases linearly (with a small slope) 

FPA Command Time Responses - Impact on Engine Location. Figure 35 shows the response to 
a step FPA command (with and without a command pre-filter) for three FPA control system designs - 
each corresponding to a distinct engine location. The three designs were implemented to achieve a single 
target FPA time response; i.e. Figure 35 shows that a similar performance profile was achievable for the 
three distinct engine locations. Figure 36 contains the corresponding elevator responses. The figure shows 
that as the engine is moved rearward (and the vehicle becomes more unstable), more elevator is required. 
As the engine is moved rearward, so does the vehicle CG. This reduces the elevator’s moment arm and it 
effectiveness - hence the need for more elevator as the engine is moved rearward. 


FPA response (w/pre-filter) FPA response (w/o pre-filter) 



Figure 35. FPA Response to Step FPA Reference Commands: Control Designs for Three Engine Locations 


37 of 45 


American Institute of Aeronautics and Astronautics 




Time (sec) Time (sec) 


Figure 36. Elevator Response to Step FPA Reference Commands: Control Designs for Three Engine Locations 


From the unfiltered elevator response in Figure 36 (right), one observes that the as the vehicle becomes 
more unstable the flexible mode effects become more pronounced; i.e. the control system cannot simulta- 
neously provide the larger bandwidth required to stabilize the more unstable vehicle while not exciting the 
flexible dynamics. The following rough rules of thumb are observed: 

• Minimum BW @ Controls: 2 times unstable pole 

• Maximum BW @ Controls: 1/2-2/3 the frequency of the first flexible mode (21 rad/s) 

The results show that a vehicle with a RHP pole between 5. 5-7. 5 becomes difficult to control with the chosen 
control system architecture; i.e. a vehicle that is too unstable may require some combination of the following: 

• a more complex control system architecture (e.g. additional notches, MIMO) 

• sensors capable of feeding back flexible mode information 

• a more rigid vehicle; i.e. larger flexible mode frequencies 


VI. B. Nonlinear Simulations 


Within this section, we consider a constant q = 2076 psf flight profile followed by a pull-up maneuver. 

The complete guidance and flight control systems may be visualized as shown below in Figure 37. To 
track constant dynamic pressure commands, a PI feedback structure (with roll off) was used for guidance: 


A guidance (s) 


k(s + z) 

Z 1 

s 

_S + Zj_ 


(41) 


where k = 1/7846, z = 15.7. z 7 is determined from the outer-loop FPA controller - it is roughly selected to 
be the FPA outer-loop’s PI controller’s zero or larger. The second term in brackets is the FPA command 
pre-filter. 


• Fuel depletion is modeled via linear decreasing mass from 6,154.1 lbs/ft to 3,049.0 lbs/ft 

• A simple PI guidance system (which processes dynamic pressure error = Q 2076 psf — (factual) is used 
to generates FPA guidance commands until Mach 8, 85 kft; 

• A step velocity reference command is used to take the vehicle along q = 2076 psf from Mach 5.7 at 70 
kft to Mach 8 at 85 kft. 

• Once Mach 8, 85 kft is reached, the pull-up maneuver is executed; To accomplish this, a FPA command 
of 3.5° is held for 100 seconds. 


Two vehicles are considered: 
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Figure 37. Inner-Outer Loop Architecture and Guidance System 


1. One vehicle is the gap optimized vehicle (blue) - faster vehicle, larger engine, larger elevator, smaller 
lower forebody inclination (more aerodynamic); 

2. The other vehicle is our nominal vehicle (with the new engine design) - slower vehicle, smaller engine, 
smaller elevator, larger lower forebody inclination (less aerodynamic); 

Figure 38 contains relevant trajectory data. Figure 39 contains the corresponding Mach, FPA, AOA, and 
FER responses. Figure 40 contains the corresponding elevator responses as well as altitude-Mach profiles. 
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Figure 38. Nonlinear Time Responses for Vehicle Maneuver: Trajectory Data 
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Figure 39. Nonlinear Time Responses for Vehicle Maneuver: Mach, FPA, AOA, FER 
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Figure 40. Nonlinear Time Responses for Vehicle Maneuver: Elevator, Altitude-Mach Profile 


In short, the figures show that 

• AOA < 12° with the peak occurring at the beginning of the pull-up 

• The vehicle’s climb to an altitude of ~ 135 kft 

• Dynamic pressure decreases to ~ 200 psf 

The benefits of the gap-optimized vehicle are as follows: 

• The maneuver is completed ~ 30% faster (due to increased thrust margin provided by the larger engine 
for the gap-optimized vehicle) 

• While the non gap-optimized vehicle control system is scheduled on dynamic pressure to accomplish 
the pull-up maneuver, this is not required for the gap-optimized vehicle’s control system. The gap- 
optimized vehicle thus requires a less sophisticated (less adaptive) control system. 

For a more detailed discussion on gain scheduling of the nominal design, see [41]. Within [80] the author 
also examines gain scheduling issues for the nominal vehicle model. 

VII. Summary and Directions for Future Work 

Summary. An engine redesign has been conducted based upon traditional as well as control-relevant 
metrics. A complete parametric study involving inlet capture area, diffuser area ratio, internal nozzle ratio, 
and nozzle exit area is presented to justify the engine redesign. Care was been taken to ensure that the 
engine parameters are feasible with respect to the geometry of the vehicle. 

Vehicle geometry studies were also conducted. These involved engine location and forebody compression 
ramp angle. Their impact on the vehicle’s trimmable envelope (region in altitude-Mach space), static (equi- 
librium) conditions at trim, and dynamic properties as they relate to control system design was addressed for 
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level flight. These vehicle geometry trade studies were then used to obtain optimized vehicle parameters for 
various performance metrics at a single flight condition (Mach 8, 85 kft): (1) maximizing lift-to-drag (L/D), 
(2) minimizing fuel consumption, FER, elevator, or AOA. The gap metric was also used to minimize dynamic 
characteristics across the flight envelope (trimmable region) . This was done with the intention of simplifying 
control system design. Finally, a simple (classic inner-outer loop decentralized) control system architecture 
was used to illustrate how vehicle parameter selection impacts control system design. Specifically, it was 
shown how a gap-optimized vehicle can simplify control system parameter scheduling. In short, the paper 
offers contributions to control- re levant vehicle design. 

Conclusions. This paper has shown how a nonlinear (3DOF + flexibility) longitudinal model can be used 
for analysis, conducting vehicle trade studies, performing vehicle optimization, and control system design. 
It was shown that a gap-optimized vehicle can simplify the control system design/scheduling process. 

Directions for Future Work. The work presented in this paper provides motivation for conducting 
comprehensive trade studies using higher fidelity vehicle models; i.e. 6DOF + flexibility 79 . As such, the 
work motivates the development of general 6DOF tools that adequately address control-relevant modeling, 
analysis, and design issues for hypersonic vehicles during the early vehicle conceptualization/design phases. 
One specific concern will be to assess when conclusions obtained from a 3DOF model may be misleading. 
This work is currently being pursued by the team 79 . 
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